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O . Abstract 

We investigate the Q 2 evolution of parton distributions at small x values, obtained in 
the case of flat initial conditions. The contributions of twist-two and (renormalon- 
Q ■ type) higher-twist operators of the Wilson operator product expansion are taken 

into account. The results are in excellent agreement with deep inelastic scattering 
CN| ! experimental data from HERA. 

Key-words: Quantum Chromodynamics, the deep-inelastic scattering, structure 
j> 1 function, parton distribution, twist. 

G\ ' 
On 
(N 



O 
O 



43 



1 Introduction 



The measurements of the deep-inelastic scattering structure function F2 in HERA [|T], |2[ 
have permitted the access to a very interesting kinematical range for testing the theoretical 
ideas on the behavior of quarks and gluons carrying a very low fraction of momentum of 
the proton, the so-called small x region. In this limit one expects that non-perturbative 
effects may give essential contributions. However, the reasonable agreement between 
HERA data and the next-to-leading order (NLO) approximation of perturbative QCD 
that has been observed for Q 2 > lGeV 2 (see the recent review in 0) indicates that 
perturbative QCD could describe the evolution of structure functions up to very low Q 2 
values, traditionally explained by soft processes. It is of fundamental importance to find 
out the kinematical region where the well-established perturbative QCD formalism can 
be safely applied at small x. 

The standard program to study the small x behavior of quarks and gluons is carried 
out by comparison of data with the numerical solution of the DGLAP equations by fitting 
the parameters of the x profile of partons at some initial Ql and the QCD energy scale A 
(see, for example, @, H). However, if one is interested in analyzing exclusively the small 
x region (x < 0.01), there is the alternative of doing a simpler analysis by using some of 
the existing analytical solutions of DGLAP in the small x limit (see for review). This 
was done so in Ref. Qj-H where it was pointed out that the HERA small x data can be 
interpreted in terms of the so called doubled asymptotic scaling phenomenon related to 
the asymptotic behavior of the DGLAP evolution discovered in many years ago. 

Here we illustrate results obtained recently in J7| and demonstrate some (preliminary) 
results of ||, where the contributions of higher-twist operators (i.e. twist-four ones and 
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twist-six ones) of the Wilson operator product expansion are taken into account. The 
importance of the contributions of higher-twist operators at small-x has been done in 
many studies (see WW). 



We would like to note that the results of p] are the extension to the NLO QCD 
approximation of previous leading order (LO) studies JTD|, |(|. The main ingredients are: 

1. Both, the gluon and quark singlet densities are presented in terms of two compo- 
nents ('+' and '— ') which are obtained from the analytical Q 2 dependent expressions of 
the corresponding ('+' and '— ') parton distributions moments. 

2. The '— ' component is constant at small x, whereas the '+' component grows at 
Q 2 > Ql as ~ exp (a NLO ), where 



o-nlo = 2y (d + s + D + p)lnx, 

and the LO term d + = —12/p and the NLO one D± = d±± + d±Pi/Po with d±± = 
412//(27/9o)- Here the coupling constant a s = a s /(47r), s = ln(a(Ql) / a(Q 2 )) and p = 
a (Qo) ~~ a (Q 2 )i Po and /3i are the first two coefficients of QCD /3-function and / is the 
number of active flavors. 



2 Basic formulae 

Our purpose is to show the small x asymptotic form of parton distributions in the frame- 
work of the DGLAP equation starting at some Qq with the flat function: 

fl 2 {Ql) = A a ( hereafter a = q,g), (1) 

where f^ 2 are the leading-twist parts of parton (quark and gluon) distributions multiplied 
by x and A a are unknown parameters that have to be determined from data. Through 
this work at small x we neglect the non-singlet quark component. 

We would like to note that new HERA data show a rise of F2 structure function 
at low Q 2 values (Q 2 < lGeV 2 ) when x — > (see Fig. 2, for example). The rise can be 
explained in a natural way by incorporation of higher-twist terms in our analysis (see the 
part 2.2). 

We shortly compile below the main results found in 0, [| at the LO approximation 
(the leading-twist results at the NLO approximation may be found in J7J). The full 
small x asymptotic results for parton distributions and F 2 structure function at LO of 
perturbation theory is: 

F 2 (x,Q 2 ) = e-f q (x,Q 2 ) (2) 
f a (x,Q 2 ) = f+(x,Q 2 ) + f-(x,Q 2 ), (3) 

where the '+' and '— ' components f^(x, Q 2 ) are given by the sum 

Q 2 ) = fl^Xi Q 2 ) + fa^ Q 2 ) (4) 

of the leading- twist parts /^ 2,± (x, Q 2 ) and the higher- twist parts /^ T,± (x, Q 2 ), respectively. 
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2.1 The contribution of twist-two operators 

The small x asymptotic results for PD, fa 2,± 
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A g + -A q y G {a) e~ d ^ s + 0(p) , 

{(A g + ~A^ph{a)e-~ d ^> + 0(p) 

-~A q e~ d -V s + 0(x), 
A q e- d -^ s + 0(x) , 



where d+(l) = 1 + 20 f/ (27 fa) and = 16 f/ (27 fa) are the regular parts of rf H 

d_ anomalous dimensions, respectively, in the limit n — > 1 Q. The functions I u (v = 
are related to the modified Bessel function I v and to the Bessel function J u by: 



iM = 

The variables a and p are given by 

a = 2\J \d + sln(x)\ , p 



if s > 
if s < 



(5) 
(6) 
(7) 



and 
0,1) 



(9) 



\d+s\ 



a 



\ ln(l/x) 2ln(l/x) 



(10) 



2.2 The higher-twist contributions 



Using the results in 



(which are based on calculations [12, |13[), we show the effect of 



higher-twist corrections in the renormalon case (see recent review of renormalon models in 
[|14|1 ). We present the results below making the following subtitutions in the corresponding 
twist-two results presented in Eqs.(f|)-(||): 

f^ + (x,Q 2 ) (see Eq.(§) - f^ + (x,Q 2 ) by 
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A n 



16/ f A? (2~ 



K 9a (f)-^(f 



8A|^ 
7 Q 4 
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In(a) 
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Q 2 



h(a) 



where A 2 a and A 4 . a are magnitudes of twist-four and twist-six corrections and 



K g9 (f) 

fI 2 >+(x,Q 2 ) (see Eq.(§ 
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/f<+(:r,Q 2 )by 



60 



4 F ^ 4 128//A 2 ja /2 f , 

A a ph(a) -> A a • -7^7 <; I -Ii(o-J 



45/? 2 1 Q 2 Vp" 



'A? 



Q 1 



Ia(a) 



7 Q 4 



'A 2 ' 



K ga (f) - ^ - In Q2 



^rom now on, for a quantity k(n) we use the notation k(n) for the singular part when n 
k(n) for the corresponding regular part. 
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Fig. 1 
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Figure 1: The structure function F2 as a function of x for different Q 2 bins. The exper- 
imental points are from HI The inner error bars are statistic while the outer bars 
represent statistic and systematic errors added in quadrature. The dashed and dot-dashed 
curves are obtained from fits (based on leading-twist formulae) at LO and NLO respec- 
tively with fixed Ql = 1 GeV 2 . The solid line is from the fit at NLO giving Ql = 0.55 
GeV 2 . 
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From Eqs.(|TT|)-(|i4|) one can notice that the higher-twist terms modify the flat condition 
Eq. ([[]). They lead to a rise of parton dostributions and, thus, F 2 structure function at 
low value gjj, when x — > 0. This is in agreement with HERA data 0, as it is shown in 
next section. 
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Fits to H195 data. Also shown ZEU95 data. dashes = LO solid = LO + HTren. nf = 4 




Figure 2: The structure function F% as a function of x for different Q 2 bins. The ex- 
perimental points are from HI and ZEUS @. The inner error bars are statistic while 
the outer bars represent statistic and systimatic errors added in quadrature. The solid 
curves are obtained from fits at LO, when contributions of higher-twist terms have been 
incorporated. The dashed curves are from the twist-two contributions alone. 



3 Results of the fits 

With the help of the results presented in the previous section we have analyzed F 2 HERA 
data at small x from the HI collaboration (first articles in [|T], Qj). In order to keep 
the analysis as simple as possible we have fixed the number of active flavors f—4 and 
Aj^jg( n / = 4) = 250 MeV, which is a reasonable value extracted from the traditional 
(higher x) experiments. Moreover, we put A l a = A 2ja in agreement with |I5| . The initial 



scale of the parton densities was also fixed into the fits to Qq — 1 GeV 2 , although later 
it was released to study the sensitivity of the fit to the variation of this parameter. The 
analyzed data region was restricted to x < 0.01 to remain within the kinematical range 
where our results are accurate. 

Fig. 1 shows F 2 calculated from the fit (based only on leading-twist formulae) with 
Q 2 > 1 GeV 2 in comparison with 1994 HI data (first article in (l[|). Only the lower Q 2 
bins are shown. One can observe that the NLO result (dot-dashed line) lies closer to the 
data than the LO curve (dashed line). The lack of agreement between data and lines 
observed at the lowest x and Q 2 bins suggests that the flat behavior should occur at Q 2 
lower than 1 GeV 2 . In order to study this point we have done the analysis considering Qq 
as a free parameter. Comparing the results of the fits (see |7]]) one can notice the better 
agreement with the experiment at fitted Ql = 0.55 GeV 2 (solid curve) is apparent at the 
lowest kinematical bins. 

Fig. 2 shows F 2 calculated from the fit at LO (based on leading-twist and higher-twist 
formulae) in comparison with 1995 HI and ZEUS data 0]. One can observe that these 
results (solid line) lies closer to the data than the twist-two results (dashed line). We 
have done the analysis considering Qq as a free parameter. Comparing the results of the 
fits (see [[|) one can notice the better agreement with the experiment at fitted Ql = 0.61 
GeV 2 , which is close to Ql in the analysis of 1994 HI data (see Fig. 1). 
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4 Conclusions 



We have shown that the results developed recently in || have quite simple form and 
reproduce many properties of parton distributions at small x, that have been known from 
global fits. 

We found very good agreement between our approach based on QCD and HERA 
data, as it has been observed earlier with other approaches (see the review 0). The 
(renormalon-type) higher-twist terms lead to the natural explanation of the rise of F 2 
structure function at low values of Q 2 and x. The rise has been discovered in recent 
HERA experiments 0. 

As next step of our investigations, we plan to study contributions of higher-twist 
operators to relations between parton distributions and deep inelastic structure functions, 
observed, for example, in [lf| |P7 |. 
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